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We derive an extremal equation for optimal completely- 
positive map which most closely approximates a given trans- 
formation between pure quantum states. Moreover, we also 
obtain an upper bound on the maximal mean fidelity that 
can be attained by the optimal approximate transformation. 
The developed formalism is applied to universal-NOT gate, 
quantum cloning machines, quantum entanglers, and qubit 
^-shifter. 

PACS number(s): 03.67.-a 



I. INTRODUCTION 

It is well known that certain transformations are for- 
bidden in quantum theory. For example, it is impossible 
to exactly clone an unknown quantum state |jj or to pre- 
pare a qubit orthogonal to a given unknown qubit 
These transformations can be implemented only approx- 
imately ||-^]]. A convenient measure of the quality of 
such an approximate transformation is the mean fidelity 
and the transformation is optimal if it reaches the maxi- 
mum possible fidelity. 

In quantum theory, the class of allowed transfor- 
mations consists of trace-preserving completely-positive 
maps (CP-maps) £ between input and output density 
matrices. In this paper we derive a generic extremal 
equation for the optimal CP-map which approximates 
a given desired transformation. This equation may eas- 
ily be solved numerically by means of repeated iterations 
which yields the optimal CP-map. We also derive an up- 
per bound on the maximum attainable fidelity. If we find 
a CP-map which saturates this upper bound then such 
a transformation is, by definition, the optimal one. We 
illustrate the application of our formalism on universal- 
NOT gate, quantum cloning machine, quantum entan- 
glers, and quantum ?9-shifter. 

We shall consider only the unconditional transforma- 
tions which always provide an output. This should be 
contrasted with conditional transformations whose out- 
put is either accepted or rejected in dependence on the 
result of some measurement on the ancilla. Note that the 
conditional transformations can achieve higher fidelity 
than the unconditional ones, but only at the expense of 
possibly very low probability of success. 

The paper is organized as follows. In Sec. II we intro- 
duce a useful parametrization of the CP-maps and de- 
rive the extremal equation for the optimal CP-map. We 
also obtain an upper bound on the maximum fidelity. In 
Sec. Ill we briefly describe how any CP-map can be im- 
plemented as a unitary transformation on larger Hilbcrt 



space. In Sec. IV we present examples of application of 
our method. Finally, Sec. V contains conclusions. 



II. EXTREMAL EQUATION AND BOUND ON 
FIDELITY 

Suppose we would like to implement a transformation 
between pure states \ipin) £ W and |t/w) € K, 



\i>in(x)) -> IV'out(x)). 



(1) 



Here x denotes a set of numbers parametrizing all pure 
states in Hilbert space of input states H. Note that the 
dimensions of input and output Hilbert spaces may differ 
in general, dim7i ^ dim/C. It may happen that the de- 
sired transformation (|l|) cannot be implemented exactly 
because (|l|) is not a linear CP-map and we thus want 
to find a CP-map which most closely approximates the 
transformation ([j]). 

Let us begin by introducing a convenient representa- 
tion of the CP-maps. We say that map p — > £(p) is 
positive if it preserves the positivity of operators p. The 
map £ is completely positive if and only if the exten- 
sion £u®T-H' is a positive map for any Hilbert space Ti', 
where X is an identity map. Any CP-map can be repre- 
sented by a positive operator Q . Consider the maximally 
entangled state on Ti® 2 , 



dimW 

\f)= E 

3=1 

and define the operator 

X = £h®Zh{\<p){<p\)- 



(2) 



(3) 



Note that x acts on a Hilbert space 7i®K.. It is easy to 
show that the CP-map p ou t = £(pm) can be written as 



Pout = TrnixpJn ® lx:], 



(4) 



where T stands for the transposition and 1*; denotes the 
identity operator on JC. The requirement that the CP- 
map £ should preserve the trace imposes the following 
constraint on x, 



Tr/c[x] = 1«. 



(5) 



We would like to quantify how well the CP map x 
approximates the desired transformation (|l]) . To this end 
we define the mean fidelity as 



1 



F= J dx(ip out (x)\£ (\t m {x))(^ in {x)\)\^out(x)), (6) 

where dx denotes the proper measure on space of pure 
states \ip- m {x)). 

With the help of Eq. (^) we may rewrite the expression 
(El as 



F = Tr[ X R], 



(7) 



where the positive operator R acting on TL K. is given 
by 



R 



doS(\ip ia (x))(ip in (x)\) ® |Vout(£)X^out@)|. (8) 



Our task is to find a trace-preserving CP-map x which 
maximizes the fidelity ([?]) . We take into account the con- 
straint ([5J) by introducing an operator Lagrange multi- 
plier A = A ® 1/c and we look for the maximum of the 
functional 

F[ x ] = Tr\xR}-Tr\xM- (9) 
We expand \ in eigenstate basis 

3 

and rewrite the functional (^) as 

3 

A variation of F[x] with respect to (iTj\ yields the ex- 
tremal equations, 



(R-A)r j \n j )=0. 



(11) 



We note that recently a similar equation has been de- 
rived for elements of positive operator valued measure 
representing an optimal quantum measurement fl!o|fLl[ . 

We multiply Eq. (|ll|) by (nj | and sum over j. After 
some manipulations we obtain 



X = A 



(12) 



Further we take Hcrmitian conjugate of this formula, x — 
\ /i'A ' . and insert it back to the right-hand side of Eq. 
(O) . Thus we arrive at a symmetrized extremal equation 



X 



(13) 



The Lagrange multiplier A = A ® lie can be determined 
from the constraint (||) which provides expression for A, 



A = (Tr K [R x R] 



1/2 



(14) 



We fix the square root by postulating that A is positive 
Hermitian operator. The system of coupled nonlinear ex- 
tremal equations (O) and dTJI) can be conveniently solved 



by means of repeated iterations, starting from some ini- 
tial 'unbiased' CP-map, for example a map which trans- 
forms every input density matrix to the maximally mixed 
state on /C, p- m — > l;c/dim/C. Note that the iterations pre- 
serve the positivity of x an d the constraint (|^) is exactly 
satisfied at each iteration step. 

From the formula (0) we can obtain an upper bound 
on the optimal fidelity F, 



F < Tr[ X ]i? max = dimHi?„ 



(15) 



where i? max is the largest eigenvalue of the operator R. If 
we find a CP-map which reaches the upper bound on fi- 
delity ( |l5| ) then such transformation is, by definition, op- 
timal one. We shall use this theorem below to prove the 
optimality of quantum cloning machine and universal- 
NOT gate. 



III. IMPLEMENTATION OF CP-MAPS 

Before turning to explicit examples of application we 
should briefly comment on the possibility of experimental 
realization of the CP-map x- If holds that any trace- 
preserving CP-map can be accomplished as a unitary 
evolution on an extended Hilbert space. Any CP-map 
can be written in the form of Kraus decomposition p2| 



Pout 



JT 



and the trace-preservation condition gives 



£44 



i 



H ■ 



(16) 



(17) 



Rewritten in terms of the matrix elements A^J = (k\Ai \i) 
this constraint reads 



E 

k,l 



(18) 



The number of necessary operators Ai is equal to the 
number of nonzero eigenvalues of matrix x which we de- 
note by C. In fact, the operators Ai may be associated 
with the eigenstates \tti) of the operator x, 



A (i) _ 



ri(k\(i\n) 



(19) 



where \i) G TL and \k) £ K, are states in input and output 
Hilbert spaces, respectively. 

In order to implement x as a unitary transforma- 
tion, we must work in Hilbert space with dimension 
D = Cdim/C, where we define the transformation 



010) 



k,l 



\k)\l). 



(20) 



Here |0) denotes an initial state of the ancilla and |Z) are 
orthogonal output states of the ancilla. With the help of 
the formula (|l8|) it is easy to check that the transforma- 
tion pfl) is indeed unitary. 
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IV. APPLICATIONS 

Having established the general formalism, we may pro- 
ceed to explicit examples. In all the examples, H is a 
Hilbert space of single qubit, all pure states in H can 
be visualized as points on the surface of Bloch sphere 
parametrized by angles i? and 0, 



JV 



|^(^)>=cos^|0> + e ^sm^|l>, 
and the proper integral measure is 



dx 



1 

4?r Jo 



sin i9 d-d d(j). 



(21) 



(22) 



A. Universal NOT gate 

Suppose that we have N copies of an unknown qubit 
\ip) a and we would like to prepare an inverted qubit \ip±), 
(t/j±\ip) — 0, and 



|^)=sin||0)-e^cos||l) 



(23) 



Let us denote the input state consisting of N identically 
prepared states \ip) as 



= |</>>i...|V>> 



(24) 



Since it is impossible to exactly invert an unknown qubit, 
the mapping 



(25) 



can be carried out only approximately H. The input 
Hilbert space is a bosonic subspace of H.® , i.e. a Hilbert 
space 7if N spanned by totally symmetric states of N 
qubits. Let us denote by \N, k), k = 0, . . . , N, a totally 
symmetric state where k qubits are in state |0) and N — k 
qubits are in state |1). Thus we can write 



^W-^ (cos^ 



fc=0 

We calculate the operator .Rnot as 



l) 



N-k 



\N, k). 
(26) 



Rnot = — wntidd d0(|i/>}M) ® |^x)<^x| 
47r Jo Jo 

After a straightforward integration we arrive at 



^NOT - £ ,„ N 1 ^t 1 o^ ! ^ ■ ) ( V °> < ( 1 : 



k=0 



(7V+ 1)(/V + 2) 



N 

£ 

k=0 
N 



k + 1 



(JV + 1)(jV + 2) 



\N,k)(N,k\ 



y k{N -k + 1) 
^ (N+l)(N + 2y A 



Y. t iN -t +1 l \^-i)(N, k \®\i m . 



(27) 



The matrix -Rnot is block diagonal and the calculation 
of the eigenvalues of -Rnot boils down to evaluation of 
eigenvalues of 2 x 2 matrices. The largest eigenvalue reads 
R max = l/(JV+2) and it is N + 2-fold degenerate. Taking 
into account that diiriH;?^ = TV + 1 we get immediately 
fromEq. (§J) the upper bound on the fidelity of U-NOT 
gate, 



-Fnot < 



N+ 1 
iV + 2' 



(28) 



Since the fidelity of U-NOT gate proposed by Buzek et 
al. reaches the bound (JV + l)/(N + 2) it is optimal [[§. 
Note also that our determination of the upper bound on 
fidelity of the U-NOT gate is conceptually similar to the 
derivation of the fidelity of the optimal quantum mea- 
surement @|3[. 

When we insert the operator ( p7j ) into extremal Eqs. 
( |l3| ) , (jl4|) , we may numerically calculate the optimal CP- 
map. To be explicit, let us consider the simplest case of a 
single input copy, N = 1. The input and output Hilbert 
spaces are now identical and the corresponding operator 
J?not can be written as 

%ot = ~(|01><01| + |10)(10| + |*-)(*-|), (29) 



where 



|*_> = -^(|00>-|11». 



(30) 



Notice that -Rnot is proportional to the projector to 
the subspace spanned by vectors |01), |10), and 1$-). 
Obviously, the eigenvalues of -Rnot are 1/3 and and 
-Fnot < 2/3. 

Due to the specific structure of the operator -Rnot, 
there exist infinitely many solutions to extremal equa- 
tions ( |l3| ) and (|TJ). In the numerical calculations, the 
result depends on the initial choice of X- I n & H cases, 
we obtain an optimal CP-map which reaches the bound 
2/3. If we start to iterate from the multiple of identity 
operator, \ — 1/2, then after a single iteration step we 
obtain the universal-NOT gate whose fidelity is state- 
independent m, 
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2pn + poo -Poi 
-pio 2,000 + pii 



(31) 



B. Quantum cloning machine 

An ideal 1 —* N cloning machine would prepare TV 
exact clones of an unknown qubit \ip), 



(32) 



We restrict ourselves to the symmetric quantum doners 
which produce N identical approximate copies. In this 
case the output Hilbert state is the bosonic space H® N 

and the ideal output state l^) is given by Eq. (ffit). The 
construction of the operator R c is straightforward, 



1 



2jt 



which yields 

N 



k+i 



k=0 



(JV + l)(JV + 2) 



N 

£ 

fc=0 
JV 



N-k + 1 



+£ 



{N + l)(N 

y/k{N-k + l) 

< (N + l)(N + 2) 



\0)(0\®\N, k)(N,k\ 



-|1><1|®|W, k)(N,k\ 



|0)<l|®|iV, k)(N,k-l\ 



+£ 



y/k(N-k + l) 

^ (N+l)(N + 2) 



\l){0\®\N,k-l)(N,k\. 



(33) 



The largest eigenvalue of the block-diagonal matrix R c 
reads 1/(N + 1)- Consequently, we have the constraint 
on maximum cloning fidelity, 

9 

F c < 



N + l 



The symmetric 1 — > N cloning machines proposed by 
Buzek and Hillery H and Gisin and Massar Q] saturate 
the bound F c = 2/(A r +l) and are thus optimal, as shown 
earlier in Refs. |MpH8j . 

For N = 2 we solved numerically the extremal Eqs. 
( |l3| ) and ( |l4| ) by means of repeated iterations and we 
obtained a CP-map which describes the optimal univer- 
sal quantum doner designed by Buzek and Hillery ||. 
The convergence was very fast and we obtained the de- 
sired CP-map with 10-digit precision after several tens 
of iterations. Since the optimal CP-map is unique the 
numerical solution does not depend on the choice of ini- 
tial CP-map xo j provided that all eigenvalues of xo are 
positive. 



C. Quantum entangler 

Recently, Buzek and Hillery investigated the pos- 
sibility of entangling an unknown qubit with a qubit 
in known state |0). They showed that the following en- 
tangling transformation 



|V0^W}|o} + |oM) = IVw}, 



(34) 



cannot be carried out exactly. The ideal output state 
|Vw) can be expressed as 



IVw(tf») 



and 



V2cos(iV2)|00) +e ^,^i sin(lV2)| , I'- 
v /l-|-cos 2 (zV2) 



|* + ) = — (|01> + |10)). 



(35) 



(36) 



Buzek and Hillery constructed a universal entangler 
whose state-independent fidelity is F « 0.946 [[uj. Here 
we show that this entangler is not optimal and we find 
an entangling machine which achieves a higher fidelity. 

We proceed along the same lines as in the two previ- 
ous examples and determine the operator i? e nt- After a 
straightforward integration we find, 

4„t = (21n2-l)|0)(0|®|00)(00| 
+(3-4m2)|l)(l|<g>|00) (00| 
+ (3/2- 21n2)|0)(0| ® |*+)(*+| 
+(41n2-5/2)|l)(l|® |* + )<* + | 
+V2(3/2- 21n2)|0)(l| ® |00)(*+| 
+V2(3/2- 21n2)|l)(0| ® |*+)(00|. (37) 

Since the largest eigenvalue of R cnt is 1/2 we have from 
( |l5| ) that F on t < 1- Although the upper bound 1 cannot 
be reached we can get closer than in Ref. [Q. When 
we numerically solve the extremal equations for optimal 
CP-map, we find that the optimal map is a unitary trans- 
formation on the space of two qubits, where the second 
qubit is initially prepared in the state |0), 



I* 



|00) - |00), 1 10) 
For an arbitrary input qubit ( ^l| ) we thus have 

| ^> |0) ^ cos 1 100> + sin 1 . 
The fidelity of this transformation depends on 
[V2cos 2 (iV2) + sin 2 (tf/2)] 2 



FentW = 



1 + cos 2 (tf/2) 



(38) 



(39) 



(40) 



We plot the fidelity F cn t(i?) in Fig. 1. The mean fidelity 
reads 
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U 

FIG. 1. Fidelity F e nt($) of the entangling machine. 



iterate Eqs. ( |l3| ) and (|l4|), we obtain the optimal CP- 
map which reaches the maximum fidelity 1/3. It turns 
out that the optimal transformation is strikingly simple. 
For any input state, the entangler should prepare the 
same mixed output state 

Pout - I (|00)<00| + |* + X*+| + |H)(11|) . (46) 

The fidelity is 1/3 and does not depend on the input 
state. Moreover, it can be easily shown with the help of 
the Peres-Horodecki criterion |l5|]l6| ] that the state jiq ) is 
separable. Thus we conclude that this optimal universal 
"entangler" exhibits rather poor performance because it 
does not prepare an entangled state. 



F ont = 3a/2 - 7/2 + (6 - Ay/2) ln(2) w 0.9805, (41) 

which is clearly higher than the fidelity 0.946 obtained 
in 0. Since F ont (i9) is state dependent, the entangler is 
not universal. Nevertheless, for any state the fidelity 
of the transformation fl39| ) is higher than the fidelity of 
the universal entangler proposed in E3 and the minimum 
of the fidelity ©is 



F min = iV2(V2 - 1) 



0.9706. 



(42) 



It is easy to see why the non-universal entangler out- 
performs the universal one. Since we a-priori know the 
second state |0), this information breaks down the sym- 
metry, and there is no reason to expect that the uni- 
versal entangler should be the best one. Indeed, there 
exists a preferred basis spanned by |0) and its orthogonal 
counterpart |1). Note that the transformation (|3|) real- 
izes ideally the desired transformation (63) for the basis 
states |0) and |1). 

Let us also study another type of quantum entangler 
considered in jl3]. Suppose we have an unknown qubit 
\tp) and we would like to prepare a maximally entangled 
state 



^(IV>>|V>±) 



IV>x)|V». 



(43) 



On inserting the desired output state ( (43| ) into (||), we 
find 

R± = (i 2 ®i 3 + ^ V a 2 j9a a j), (44) 

o 1 — ' 



where dj denote Pauli matrices and the subscripts 1 and 
2, 3, refer to input Hiibert space Ti. and output space 
K = TC®TC, respectively. The largest eigenvalue of R±, 
which is equal to 1/6, limits the fidelity of the entangler 



F ± < 1/3. 



(45) 



This result was already obtained in |Q but our proof 
is much simpler than that presented in 114]. When we 



D. Optimal ^-shifter 



Recently, Hardy and Song |17|,|1 
proximate implementation of single-qu 



studied an ap- 
Dit transformation 



where the angle & is shifted by a € [0, 7r], 
1^(0, </>)) - |#? + a,0)>. 



(47) 



Hardy and Song found that the optimal universal CP- 
map whose fidelity does not depend on •d is identity map 
for a < 7r/2 and a U-NOT gate for a > ir/2. Remarkably, 
the lowest fidelity F = 1/2 is achieved for a = ir/2. 
However, the universal transformation is not the best 
one and one can find a CP-map which for certain range 
of angles a approximates the transformation (|47| ) much 
better (i.e. attains much higher mean fidelity) than the 
optimal universal transformation |I§| ] . 

In what follows we show that our approach naturally 
and straightforwardly leads to the optimal non-universal 
CP-map. On inserting the output state ( [i"7j ) into Eq. (j^) 
we obtain the operator R$(a), 

Ma) = Q + -Lcosa-^sinaj |00)(00|, 



1 

12 
1 

12 
1 

12 



cos a 



cos a 



16 



sina |01)(01|, 



16 



sina |10)(10|, 



16 



sina |11)(H|, 



— cosa(|00)(ll| + 111) <00|). 
6 



(48) 



We have solved numerically the extremal Eq. (|T^), jl^ ) 
for various values of a. An analysis of the structure of 
resulting CP-maps reveals that they represent a simple 
damping process and we can make for them the following 
ansatz, 
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-> cos 2 /3|0)(0| +sin 2 /3|l)(l|, 
|1>(0| ->cos0|l)(O|, 
|O)<l|->cos0|O)(l|, 

|1><1| - |1><1|. (49) 
The mean fidelity is a function of a and (3, 

1 cos a , o 7T , 
F(a, 0) = - H — (cos 2 (3 + 2 cos f3) + - sin a sin 2 /3. 

2 6 8 

(50) 

The optimal (3 is obtained from the extremal Eq. 
i) 



dp 



F(a,0)=O. 



The solution 

exists for all a. A second solution 



, 37T 

cos p2 = [ — tan a — 1 



exists only if 



a > arctan — = a . 

3lT 



(51) 



(52) 



(53) 



(54) 



The second root always leads to higher mean fidelity than 
the first one. Thus for a < ao it seems optimal to apply 
an identity operation. However, as soon as the angle shift 
a overcomes the threshold ao one should rather apply a 
damping process (49) with angle (3 given by Eq. (53). 
The resulting fidelity F(a) can be expressed as 



F(a) 



|(1 + cos a), a < ao, 
F(a,(3 2 ), a>a , 



(55) 



and is plotted in Fig. 2. Notice that F(a) is not a mono- 
tonic function of a |l^] . In Fig. 2 we also show the upper 
bound on fidelity obtained from the maximum eigenvalue 
of R$ (a) . We can see that the fidelity ( |55| ) of the approx- 
imate ^-shifter ( fl9| ) is lower than the upper bound. This 
bound is reached only at three discrete points. The most 
trivial case is a = where the identity transformation 
achieves fidelity 1. For a = tt the U-NOT gate reaches 
maximum possible fidelity 2/3. The third point a = tt/2, 
where the fidelity of the i9-shifter reads 



(I) 



is perhaps most interesting. Since for a = tt/2 we have 
from Eq. (||) (3 2 = tt/2, the output of the tf-shifter @ 
is the state |1) regardless of the input state. This means 
that the transformation d — > ■& + tt/2 is implemented 
exactly for all pure states on the equator of Bloch sphere 
because all these states should be mapped to the south 
pole of this sphere, i.e. to the state |1). 




FIG. 2. Fidelity F(a) of the i9-shifter (solid line). The 
dashed line shows the upper bound on the fidelity obtained 
from the largest eigenvalue of R$(a). 



V. CONCLUSIONS 

We have derived an extremal equation for optimal com- 
pletely positive map \ which most closely approximates 
certain prescribed transformation between pure states. 
The nonlinear extremal equation can be conveniently nu- 
merically solved by means of repeated iterations, which 
provides very straightforward and simple means for deriv- 
ing the optimal CP-map. Moreover, we have obtained an 
upper bound on the mean fidelity which, in some cases, 
allows us to simply prove the optimality of the calcu- 
lated CP-map. We have applied our results to universal 
NOT gate, quantum cloning machine, quantum entan- 
glers, and qubit i9-shifter. The examples discussed in 
this paper illustrate that the present method is very gen- 
eral and may find further applications in quantum infor- 
mation processing, for example in design of optimal CP 
maps for quantum teleportation p9[. 
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